Fluctuation-Exchange Study of Antiferromagnetism in Electron-Doped Cuprate 

Superconductors with Disorder 
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On the basis of the Hubbard model, we extend the fluctuation-exchange (FLEX) approach to 
investigating the properties of antiferromagnetic (AF) phase in electron-doped cuprate supercon- 
ductors. Furthermore, by incorporating the effect of scatterings due to the disordered dopant-atoms 
into the FLEX formalism, our numerical results show that the antiferromagnetic transition temper- 
ature, the onset temperature of pseudogap due to spin fluctuations, the spectral density of the single 
particle near the Fermi surface, and the staggered magnetization in the AF phase as a function of 
electron doping can consistently account for the experimental measurements. 
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The cuprate high-temperature superconductors are the 
typical quasi-two-dimensional (Q2D) strongly-correlated 
electron systems. At low temperature, the electron- 
doped cuprates (EDC) such as Nd2-i,Ce:j;Cu04 (NCCO) 
and Pr2-2;Cea:Cu04 (PCCO) are in the antiferromag- 
netic ( AF) state within a wide doping range Xi 2.: J.j _4] • 
The Fermi surface (FS) evolution with electron dop- 
ing observed by the angle-resolved photoemission spec- 
troscopy experiment j5|, (SI reveals the microscopic infor- 
mation about the antiferromagnetism. For understand- 
ing the physics of FS evolution in EDC, much of the 
works have been carried out by using either the Hub- 
bard or t - J model 0, H H EE El El El • But there 
are difficulties with the existing theories in describing the 
AF phase transition and the FS evolution in a consistent 
manner. To explain the FS evolution and the staggered 
magnetization as a function of doping, most of the the- 
ories have to assume a doping-dependent Hubbard U . 
Since U is the on-site Coulomb interaction, it is hard 
to understand how it could be screened by conduction 
electrons. Though the mean-filed theory (MFT) yields 
qualitative explanation for FS evolution [3, [l^ llol > ^^^ 
predicted AF transition temperature T/v is well known 
to be too high 10, 13] . In order to reduce the magnitude 
of T/v, the self-consistent random-phase approximation 
which takes into account the effect of spin fluctuations 
has been employed, but the obtained T/v [lOj is a non- 
monotonic function of doping concentration x, contrary 
to the experimental observations. Therefore, a quanti- 
tative understanding of the AF phase transition in EDC 
remains to be a challenging issue which so far has not 
been properly addressed. 

In this work, we extend the fluctuation-exchange 
(FLEX) approach based on the Hubbard model [lj|, 
which was designed for studying the spin fluctuations in 
a two-dimensional strongly-correlated electron system, to 
investigating the antiferromagnetism below T/v in NCCO 
or PCCO. Furthermore, by incorporating the effect of 
scatterings due to the disordered dopant-atoms Ce into 



this formalism, we are able to show that the numerically 
obtained T/v, the onset temperature of pseudogap due to 
spin fluctuations, the spectral density of the single par- 
ticle near FS, and the staggered magnetization m in the 
AF phase as a function of x are all in good agreement 
with the relevant experimental measurements. Here U is 
doping independent and each dopant Ce atom is expected 
to release one electron in the CuO planes. With increas- 
ing doping, the scattering effect due to dopant atoms 
becomes more pronounced which is expected to expedite 
the suppression of the AF order. 

For describing the electron system, we use the Hubbard 
model including the impurity fleld: 
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where ^j^ — e^: — /i with e^ as the tight-binding en- 
ergy, /i the chemical potential, cl {ct ) is the electron 
creation (annihilation) operator of spin-cr(= ±1 for up 
and down spins, respectively) in state of momentum fc. 
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ct , ~ , and N is the total number of the 



lattice sites. Vi in "the last term is the impurity potential 
acting on the electron with density rii at position i. In 
NCCO or PCCO, the concentration Ni of the dopant Ce 
atoms is the same as the electron doping concentration 
X. The electron system under consideration is the Q2D 
one with weak interlayer coupling. 

We begin to work with the Green's function G. The 
formalism for determining G is to flnd out the self-energy 
S as a functional of G. In the normal state, follow- 
ing the FLEX approximation [lj|, we can easily get the 
self-energy functional in the presence of impurities. S 
is diagrammatically shown in Fig. 1. The diagrams in 
Fig. 1(a) come from the Hartree term. The first one 
includes the diffuson correction due to the particle-hole 
propagator under the impurity scattering. The second 
one contains the Cooperon correction. Figures 1(b) and 
1(c) are the FLEX diagrams with the decoration of the 
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FIG. 1: Diagrams for the self-energy in normal state. The 
wavy line represents the on-site Hubbard U. The dashed 
single line with a cross denotes the impurity scattering, (a) 
Lowest order (Hartree) terms, where a double counting (d.c.) 
without the impurity scattering should be subtracted, (b) 
ring diagrams representing the longitudinal spin and den- 
sity fluctuations, (c) diagrams for transverse spin fluctuations, 
(d) vertex corrections due to the impurity scattering for the 
particle-hole and (e) particle-particle channels, (f) self-energy 
due to the impurity scattering. 



impurity scattering. The ring diagrams in Fig. 1(b) rep- 
resent the contributions from the longitudinal spin and 
density fluctuations. The transverse spin fluctuations are 
taken into account in Fig. 1(c). The impurity insertions 
in the above diagrams are illustrated in Figs. 1(d) and 
1(e). Such insertions ensure the spherical symmetry for 
the longitudinal and transverse spin susceptibilities in the 
presence of impurity scatterings in the normal state. The 
process of the single electron scattering off impurities is 
shown in Fig. 1(f). Since the antiferromagnetism below 
T/v is also under consideration, instead of writing down 
the formula for Fig. 1, we need to extend the normal- 
state expressions to include the AF order below T^. 

The AF state is characterized by the staggered magne- 
tization m with commensurate wavevector Q. In the AF 
state, the operator n^ — > —aNm/2 is a macroscopic 

quantity. The term of q = Q in the Hubbard interaction 
in Eq. Q should be singled out, and be treated as 
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where A — Um/2 is the order parameter. Though Eq. 
(0) seems like the MFT, we here also take into account 
simultaneously the fluctuation terms q ^ Q m. the inter- 
action. Because of Eq. (O, we need to use the Green's 
function G = —{Tripj:^{T)^t (r')) defined as a 2 x 2 ma- 
trix with r/'i = (ct ,acl - ). Since the potential (± A) 

^ k<T ^ kcr' k+Qa' ^ ^ ' 

in Eq. |2Jl for opposite spin electrons has opposite sign. 
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FIG. 2; Off-diagonal self-energy E12. 



G so defined is independent of the spin. The order pa- 
rameter A can then be written as 
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where /3 is the inverse of the temperature T, and k = 
{k, iujn) with uj„ the electronic Matsubara frequency. 
From Eq. |2Jl with the impurity infiuence, the off- 
diagonal self-energy Si2(fc) is obtained as in Fig. 2. 

To extend the normal-state self-energy shown in Fig. 
1 to that of the AF state, we note that the transverse 
spin fiuctuation in Fig. 1(c) should include the collective 
modes, especially the Goldstone mode (GM). The latter 
is consistent with the gap equation ||2Jl. This consistency 
gives rise to a constraint on the form of E12, and one 
can show that Fig. 2 is the only graph for S12 within 
the FLEX approximation. In addition, the insertion of 
impurity scatterings in the transverse spin fluctuation is 
between the opposite spin electrons, while in Fig. 2 it is 
between the same spin electrons. In order to ensure the 
consistency between GM and Eq. Q, these insertions 
should be the same; so the processes due to G12 should be 
excluded. Keeping these points in mind, the expressions 
for the self-energy in the AF state can be obtained. 

To give an explicit expression for the self-energy, we 
start with the vertex corrections shown in Figs. 1(d) and 
1(e). They can be expressed by the same formula. 
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withX{q,iuJn) = ;i- ^^^ Gii(fc, iw„)G'ii(fc-t-g, iw„-f zfi^), 
where q = {q,ifljn) with fl„i the bosonic Matsubara fre- 
quency. With the vertex correction, the logitudinal Xz 
and transverse %* spin susceptibilities are given by 

Xz{q) = -s^^{q,i(^n)X+{q,iuJn), (5) 
XiM) ^ -^^'^{q,-i^n)X-{q,iuJn), (6) 

n 

xLC?) ^ « X! ^(^' i^n)Xo{q, iuJn)T{q + Q, zw„), (7) 

n 

with 

X±{q,tuJn) = 1 ^[Gii(fc)Gii(fc + q)± Gi2(fc)Gi2(fc + q)] 

k 

Xo{q,iuj,,) = 1 5][Gii(fc)Gi2(fc + q)- Gi2(fc)Gii(fc + q)] 



and Q means Q = {Q,0). Xiil?) is the susceptibility 
between the transverse spin fluctuations of the same mo- 
mentum q, while Xi2il) between that of the momenta q 
and q + Q. For X21 and x|2> we have xiiil) = Xuil) 
and X22(9) = Xii(9 + Q)- The expression for Sn reads, 

with Vcff{q) = V2(g) + Vii(g) and 

^11(9) - -U{[1 + Ux\q)]-' - 1 + Ux\q)}n, (10) 
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FIG. 3: Transition temperature Tn as function of doping con- 
centration X. The circles are the onset temperatures T* of the 
pseudogap. The solid circles ^l], triangles Q, squares Q, and 
the circles |ld| in the inset are the experimental data. The 
lines are for the eye. 



The first term in right-hand side of Eq. comes from 
the Hartree terms in Fig. 1(a). The contribution to the 
self-energy in the absence of impurity scattering is sub- 
tracted because of a double counting and the remaining 
constant absorbed in the chemical potential. S^ corre- 
sponds to Fig. 1(f). The last term is due to the density 
and spin fluctuations. At q = Q, Vn is divergent because 
of 1 -f UxiiiQ) = 0, which means the existance of GM 
and is consistent with Eq. (PJ. In the normal state with 
G22 = 0, if Vt = 0, Eq. ® reduces to the FLEX theory 
[13, UM ■ The off-diagonal self-energy according to Fig. 2 
can be written as 



Ei2(fc) = Ar(Q,zc^„)- 
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The other elements of S(fc) are given by E2i(fc) = Si2(fc), 
andI]22(fc) = Sii(fc-f Q). 

We here consider the interlayer coupling (IC). For 
NCCO or PCCO in which the CuO layers are staggered 
stacking, the AF spin coupling between the nearest lay- 
ers is determined by the local lattice distortions or any 
other asymmetry; the ideal stacking leads to zero spin- 
coupling. We here will not consider this complicated situ- 
ation in detail, but artificially introduce a weak IC in the 
following manner. Note that IC becomes significant only 
at the singularity in Eq. (fTU)) where 1 -I- UxiiiQ, 0) = 
with Q = (tTjTt, Qz). At (f — > Q, we should have 
1 + Uxii{q,0) ~ |c- {q — Q)P, where c is a constant vec- 
tor. The IC is refiected by the small z-component con- 
stant Cz- This constant controls the overall magnitude 
of T/v- However, the shape of the AF phase boundary 
should be determined predominantly by the intra-layer 
coupling. Noting that G and x's very weakly depend 
on the weak IC, we then approximate them as the two- 
dimensional functions by adding a (72-dependent term 
c1{qz —QzY to the denominator of the final expression of 



{[l-f Ux\q)]~^}ii in the first term of Eq. (UUJ and also 
to 1 -I- Uxz{q) (which vanishes at T/v) in the first term of 
Eq. @. From our numerical calculation, the magnitude 
of in-plane component of c is found to be near the order 
of unity. We then choose c1 = 10~^ in our calculation. 
Since the gz-integral is over a period with Qz centered, 
the final result is independent of Qz ■ 

For studying the AF properties in EDC, in our nu- 
merical calculation, the parameters in unit of t = 1 [the 
nearest-neighbor (NN) hopping] in the model are cho- 
sen as t' = -0.25 (next NN hopping), t" = 0.1 (third 
NN hopping), U — 8, and Vi = 2. The parameters t', 
t" and U, are the typical ones used in the literatures 
9]. With the recently developed algorithms for dealing 
with the summation over the Matsubara frequency and 
the Fourier transform of effective interaction with long- 
wavelength singularity [17l | , we have numerically obtained 
the self-consistent solution to the Green's function. 

In Fig. 3, the numerical result for Tn is presented 
as a function of x. The curve of u^ = corresponds to 
the FLEX calculation. For comparison, the experimen- 
tal data for NCCO are presented in the inset [ij, |j, y{ . It 
should be mentioned that the present result for Tm is one 
order of magnitude smaller than that of MET |lOl [iSJ . 
Even though the FLEX gives the overall lower T/v, its 
nonmonotonic behavior does not reflect the feature of 
the experimental data. In contrast, the calculation with 
impurity strength Vi — 2 can produce a monotonic de- 
creasing function T/v oi x, in fairly good agreement with 
the experiment if t is choosing to be 0.38 eV. This re- 
sult stems from a combination effect of spin fluctuations 
and impurity scatterings. For comparison, we also depict 
our result on the onset temperature T* of the pseudogap 
which is originated in the spin fluctuations and exhibits 
itself as a depression in the density of states around the 




FIG. 4: (Color online) Occupied spectral density around the 
chemical potential at doping concentrations x — 0.05, x — 
0.11, and X = 0.14, and at T/t = 0.003. 



chemical potential. Our calculated T* is consistent with 
the experimental data extracted from measurements on 
optical conductance [ig. Similar result on T* has also 
been obtained by Kyung et al. [ll| using a different ap- 
proach with doping dependent U. 

Shown in Fig. 4 are the occupied spectral densities 
-J{E)hnGii{k,E + i0+)/7r (with / the Fermi distribu- 
tion function) at the chemical potential integrated within 
the energy window of — 0.05t < E < 0.05i for three dop- 
ing concentrations at T/t — 0.003. In the AF phase, 
the original single band folds into a lower and an upper 
Hubbard bands. At a; = 0.05, there exists extended gap 
from (0,7r) to (7r,0). At a; = 0.11 and 0.14, gaps show 
up around the hot spots. For larger doping where the 
AF phase diminishes and the spectrum becomes a single 
curve. This result essentially illustrates the FS evolution 
with X, which has also been qualitatively explained by 
the two-band picture in the MFT 0, E3. EW- As seen 
from Fig. 4, the feature of the experimental result is well 
reproduced by the present calculation. 

Figure 5 exhibits the staggered magnetization m as 
a function of x at T = OK (obtained by extrapolation 
from the results of finite temperatures) and the com- 
parison with the experimental data [3|, ll9j . The exper- 
imental value for m{0) at Cu^+ is about 0.5 in NCCO 
and ^ 0.4 in PCCO |j|, which are comparable with the 
present result 0.44. The behavior of m obtained by the 
present calculation is in very good agreement with the 
experiments. Kusko et al. 7] using the MFT with a 
doping-dependcnt U have obtained m{x) which is a lin- 
early decreasing function of x with r7i = 0ata::«0.14 and 
m{0) ~ 0.4. This result agrees well with the experimen- 
tal measurements except near the doping region where 
m goes to zero. Though the dynamical-MFT calculation 
by Senechal et al. [9J yields a reasonable behavior for 
m{x)/m{0), to(0) ~ 0.7 is too large. 

In summary, we have extended the FLEX approach 
to the investigation of the AF phase in electron-doped 
cuprates. A self-consistent treatment of the effect of spin 
fluctuations in combination with scatterings of electrons 
due to the disordered dopant atoms lead to the suppres- 
sion of both m and T/y. The present approach provides 
numerical results which could consistently account for the 
experimental measurements on the doping dependences 
of the Fermi surface evolution and the staggered magne- 
tization m at low temperature, as well as the AF transi- 




FIG. 5: Staggered magnetization m as function of x at 
T = OK. The inset shows experimental data (solid and hol- 
low circles |3|, solid triangles [13) compared with the present 
calculation. 



tion temperature T/y and the onset temperature T* for 
pseudogap formation. 
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